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Abstract 

Humans routinely solve problems of immense computational complexity by intuitively 
forming simple, low-dimensional heuristic strategies |[T]l2l- Citizen science exploits this 
ability by presenting scientific research problems to non-experts. Gamification is an effec¬ 
tive tool for attracting citizen scientists to provide solutions to research problems. While 
citizen science games Foldit ||3l, EteRNA ||4l and EyeWire Q have been used successfully 
to study protein and RNA folding and neuron mapping, so far gamification has not been 
applied to problems in quantum physics. Does the fact that everyday experiences are based 
on classical physics hinder the use of non-expert citizen scientists in the realm of quantum 
mechanics? Here we report on Quantum Moves, an online platform gamifying optimization 
problems in quantum physics. We show that human players are able to find solutions to dif¬ 
ficult problems associated with the task of quantum computing |6]. Players succeed where 
purely numerical optimization fails, and analyses of their solutions provide insights into the 
problem of optimization of a more profound and general nature. Based on player strate¬ 
gies, we have thus developed a new, few-parameter heuristic optimization method which 
efficiently outperforms the most prominent established numerical methods. The numerical 
complexity associated with time-optimal solutions increases for shorter process durations. 
To better understand this, we have made a low-dimensional rendering of the optimization 
landscape. These studies show why traditional optimization methods fail near the quantum 
speed limit ITHUlll, and they bring promise that combined analyses of optimization land- 
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scapes and heuristic solution strategies may benefit wider classes of optimization problems 
in quantum physics and beyond. 

Quantum physics holds the potential of unpreeedented teehnologieal advanees in the realms 
of eomputing iflM and simulations ifTTll . To ensure funetionality, all quantum operations must be 
exeeuted near perfeetion, requiring highly optimized operations with fidelities above F > 0.999 
ttH. Given the high dimensionality of quantum optimization problems, one eould expeet these 
quantum optimal control problems to be impraetieally diffieult to solve. However, assuming 
full eontrollability, quantum optimization problems are benign sinee all loeal maxima are also 
global maxima llo . Tailored loeal optimizers, sueh as the gradient-based Krotov Algorithm 
(KA), solve these problems ffldl . 

Quantum eomputing operations must be exeeuted faster than typieal deeoherenee times to 
ensure funetionality. However, there is a shortest proeess duration with perfeet fidelity, denoted 
the quantum speed limit (QSL) 0, whieh imposes a fundamental limit on the proeess duration 
and henee quantum eomputation. With limited duration, the quantum optimization problem 
loses the favourable properties stated earlier and loeal optimizers are not almost guaranteed to 
eonverge. In this ease the prevalent and hitherto sueeessful paradigm is using multistart of sueh 
loeal algorithms [fT5l [HI . For a time-dependent Hamiltonians the QSL is eonventionally 
eomputed by assuming that it eoineides with the proeess duration for whieh the multistarting of 
loeal optimization fails lIT^ lSll. We demonstrate that this assumption is not neeessarily true. 

High-dimensional optimization problems are often solved by humans using simple heuristie 
strategies. Common examples are visual pattern reeognition IITtII and eatehi^ a flying ball 
subjeet to wind and air resistanee 0. Citizen seienee projeets sueh as Foldit [[3|, EyeWire [|51 
and Galaxy Zoo IfT^ employ these human skills to solve highly eomplex researeh problems via 
gamifieation. We ask whether eitizen seienee projeets ean be extended from these puzzle and 
pattern reeognition tasks to dynamie ehallenges, and whether this approaeh ean be implemented 
on quantum physios problems. In order to investigate these questions we oreated Quantum 
Move^ whieh presents quantum eomputing operations as games. 

The quantum oomputer arohiteeture we gamify employs neutral atoms trapped in optioal 
lattioes [flTTl . By driving the superfluid-Mott insulator transition, samples of hundreds of atoms 
have been trapped in optioal lattioes with an unpreeedented purity [fT^IIOll . Their regular spao- 
ing offers the possibility of ereating a sealable quantum oomputer, whieh is a ohallenge for other 
available arohiteetures [|2T1. Several proposals for the implementation of quantum eomputing 
in this system have been proposed, mainly using long range gates or oontaet interaetions IfTTll . 
Here we investigate an arohiteeture 0 where oontaet interaotion is aohieved by moving atoms 
on top of eaoh other using a so-ealled optieal tweezer Il22ll2^ . a tightly foeused off-resonant 
laser beam. Finding the optimal motion of the tweezer from one lattiee site to another is a diffi- 
oult problem when the available time is olose to the QSL, and it is this transfer problem that is 
introduoed to the players of Quantum Moves through different ehallenges. 

Here we present the results of one ehallenge ealled BringHomeWater (BHW). The aim of 
BHW is to move the optieal tweezer into a region where an atom is trapped in a fixed potential 
well, eolleet the atom and move it baek to a target area as quiekly as possible. Fast movement in- 
troduees exeitations in the state of the atom, whieh is deseribed by a quantum meehanieal wave 
funetion ^(x, t) and visualized as a sloshy liquid—see Fig. [^for the player view of BHW and 
Methods for a deseription of the game interfaee and player demographies. The exeitations must 
be stabilized before reaehing the target area to attain high fidelity of the underlying quantum 
proeess. 

^Found at www.scienceathome.org 
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BHW belongs to the class of one-dimensional quantum optimal control of individual atoms 
with one or two control parameters. This class of problems is of great interest, has been studied 
extensively within the past decade (see e.g. Refs. iHllllMl) and is highly relevant for the exper¬ 
imental platforms in Refs. [[2011251 . We believe our gamification strategy is extendable to these 
problems and other physical interactions including, but not limited to, many-body dynamics in 
optical lattices [l26l and Bose-Einstein condensates [|25l . 

Recall that the standard approach for solving a problem such as BHW is to use a multistart 
of tailored local optimization algorithms, such as the Krotov algorithm (KA) [IT4l[T6l [8l. The 
choice of an initial seed is a central challenge in such complex optimizations. In BHW the high 
dimensional optimization space can only be searched sparsely, and we used multiple methods 
for creating initial seeds for the KA. The most successful method KASS employs linear combi¬ 
nations of sinusoidal functions as seed trajectories (see Methods). Using the KA on these seeds, 
high fidelity solutions are readily found for process durations longer than T = 0.40 (units de¬ 
fined in Methods). Iteratively, solutions with shorter durations are found by contracting the 
solutions with a slightly longer duration and using them as seeds for the KA. This procedure 
is called a sweep, and it traces out entire families of solutions. The best results of KASS are 
displayed in red in Fig. |^. This method locates an ostensible QSL at = 0.29 after ap¬ 
proximately 7.4 ■ 10® trials. This KASS optimization formed our most successful bare computer 
optimization method, outperforming also the acclaimed CRAB algorithm [|27]| . 

KASS and CRAB fit the prevalent paradigm of multistarting of local optimizers, while 
global optimizers are rarely used in quantum optimal control. To investigate the BHW problem 
using a global optimizer we chose the so-called differential evolution algorithm (DE) due to its 
demonstrated success in quantum problems [|28|. For duration T = 0.40, DE performed worse 
than the Krotov optimization with identical computational resources. The poor performance 
of DE can be attributed to the very high dimensionality of the optimization space, and to the 
scarceness of good solutions found therein. 

Players were trained in a series of introductory levels before reaching the scientific chal¬ 
lenges. Training levels equip the players with a range of skills; successful solutions of the 
scientific challenges require a holistic combination of these skills. In total, all Quantum Moves 
games have been played about 500,000 times by roughly 10,000 players. BHW is the most 
played scientific level (approx. 12,000 plays by 300 players) and is therefore analysed in detail 
here. Player results span many different fidelities and process durations—see the dots in Fig. |^. 
Remarkably, players trace out a region in the vicinity of Tqg™ very similar to that of the best 
tained numerical results, despite the fact that the numerical optimization used roughly 100,000 
times more trials than the players. For short durations, players find even better solutions than 
the numerical optimization, albeit with imperfect fidelities. 

This result inspired us to introduce a powerful hybrid Computer-Human Optimization (CHOP) 
scheme, in which we use the players’ intuitive solutions as seeds for the local numerical opti¬ 
mization. CHOP was applied to the top 70 % of the player solutions with durations shorter than 
T = 0.40. The results of these player-seeded optimization sweeps are shown in Fig. ig>. It 
turns out that the CHOP solutions bunch in groups, which we denote as clans and discuss in 
detail later. In Fig. [^we illustrate the best solution families of the two dominant clans (blue 
and yellow curves).^rhe QSF found by CHOP is T^g™ = 0.20, a vast improvement of the 
value obtained from the bare computer optimization using multistarted local optimization. This 
result clearly demonstrates that BHW is a quantum control problem without the benign prop¬ 
erties discussed previously and questions the common assum ptio n that the QSF can be found 
numerically as the duration for which local optimization fails [11611 . 

To better understand the ability of the players and CHOP to give better solutions than KASS, 
we examine the clustering of the CHOP solutions. We introduce a measure describing the 


3 


distance between two solutions for a partieular proeess duration T: 

= {fjk\fjk)dt, (1) 

where \fjk) = — \^jJk{x,t)) is the differenee between two wave funetions evaluated 

along the path. Figure shows a distanee map of the CHOP solutions. Two dominant elans 
stand out as distinet regions in the distanee map—See Methods for details on the identifieation 
of elans. The solutions of the BHW transport problem eorrespondirm to these two elans are 
shown in Fig. |^. Using the first solution, marked in yellow in Fig. pp, the atom is eolleeted 
by tunneling the wave funetion into a tweezer potential plaeed on tne left hand side of the 
statie potential. In the seeond elass of shoveling solutions, marked in blue in Fig. [^, the 
tweezer is moved past the position of the atom so that the overlap of the tweezer and the statie 
potential forms a strong potential gradient aeeelerating the atom towards the target. Fast non- 
adiabatie solutions must spread the wave funetion into different energy eigenstates. In Fig. 
we see that the solutions in both elans populate many different eigenstates. The tweezer must 
then be shaken periodieally when approaehing the target area to bring the atom baek to the 
desired ground state—see Methods for more details. A priori, it was not obvious that two 
strategies, eorresponding to distinetly different physieal phenomena, should exist for this ki nd 
of a problem. It is also worth stressing that players explore both solutions despite having no or 
little prior knowledge of quantum meehanieal phenomena sueh as tunneling. 

The diffieulty of a partieular optimization problem ean often be assessed with knowledge 
of the topology or ruggedness of the so-ealled optimization or fitness landsea pe. This has been 
established for problems in quantum optimal eontrol [l29ll and in other fields 11^ . In this land¬ 
seape the quality of the solution for eaeh set of the eontrol parameters is represented as the 
height. If many global optima are distributed aeross the landseape, loeal gradient-based seareh 
is often suffieient to find the highest peak. On the other hand, if the landseape is very rugged 
and eontains many loeal maxima, loeal methods will in general fail. In order to understand our 
problem in this terminology we did a dimensional reduetion of the data from the player seeded 
solutions and the bare numerieal optimization. We assign points to individual solutions in a 
two-dimensional spaee sueh that the Euelidean distanee between two points approximates the 
Manhattan type distanee between two solutions in the full high dimensional optimization spaee 
(see Methods). The height of the landseape quantifies the fidelity of a solution. 

For proeess durations T = 0.40 and T ~ 0.17, above and below , the landseapes are vi¬ 
sualized in Fig. 1^ andj^, respeetively. Green areas in the landseape mark the player-generated 
CHOP solutions. For the longer duration (Fig. 1^), global maxima indeed spread aeross the op¬ 
timization landseape, explaining the sueeess ofKJOtov-based methods. For the shorter duration 
(Fig. 1^) all high-fidelity solutions lie in the green CHOP region. This explains the failure of the 
Krotov-based methods and diffieulty in loeating the true Tqsl, sinee the global optima are no 
longer spread aeross the landseape. Thus Fig. and|^ demonstrate the dramatie ehange in the 
landseape as the duration is deereased. It is worth stressing that the CHOP region is tiny in the 
high-dimensional landseape and therefore easily missed even by elaborate seeding strategies. 
CHOP outperforms the numerieal optimization beeause players are able to heuristieally identify 
the regions of high fidelity in the high-dimensional optimization landseape, thereby finding the 
best seeds for Krotov-type loeal algorithms. 

One long term vision of our work is to eireumvent the need for gamifieation by learning 
how players form their sueeessful low-dimensional heuristie strategies, and to ineorporate this 
into autonomous optimization algorithms. As a first step in this direetion, we introduee here 
a Heuristieally Initialized Loeal Optimizer (HILO) algorithm. HILO parametrizes the player 
solutions in a low-dimensional subspaee, whilst retaining the main features of good seeds. A 
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local search algorithm can then move beyond this subspaee to find optimal solutions. More 
specifieally, inspired by the player solutions we eonstructed a three-dimensional parameteriza¬ 
tion consisting of moving the tweezer i) right, ii) slowly left, and iii) quiekly left (see Methods 
for details). Paths from this three-dimensional spaee for T = 0.15 were used as seeds for the 
KA and iteratively applied to shorter and longer durations with a sweep (purple eurve in Fig. 
1^). Seeds were taken from the low-dimensional subspaee using a simple direet seareh. HILO 
nnds the lowest QSL at = 0.19, outperforming even the best CHOP strategies. The so¬ 
lutions from HILO are shown as a landseape in Fig. |^. Initialized in a low-dimensional seed 
spaee, HILO effieiently explores a smaller but more optimal volume of the global optimization 
spaee. Only parametrizations whieh aeeurately eapture the nature of effioient solution strategies 
at short durations will lead to effioient optimization. The intuition gained from the players was 
pivotal as the parametrization used in HILO emerged from the CHOP solutions and the physi- 
oal interpretation of the player-based solution strategies. This makes CHOP a orueial preeursor 
to HILO. We stress that any optimal oontrol effort requires a substantial amount of optimiza¬ 
tions of different seeds to find good solutions. After this initial work, HILO could be readily 
implemented by first applying our elustering methods to identify potential elans. If successful, 
analysis of the elans would allow for a formulation of the low-dimensional parametrization used 
in HILO. 

Finally, we diseuss another central eomponent of quantum optimal control theory, namely 
the dependenee of the fidelity on the process duration around the QSL. This has previously 
been associated with a universal-type sin^(|T /Tqsl) behavior dH. Here we do not observe this 
behavior (see Extended Data Fig. [^i and Methods for details). We attribute this to variations 
in the so-ealled direet Hilbert speM dH. The variations in our problem may arise from the 
sequential nature of the good strategies, sueh as the three steps in the HILO-parametrization, 
where each sequence has a different scaling of fidelity with time. 

Using the gamified interface in Quantum Moves, players with little or no training in quan¬ 
tum physios not only provided high quality solutions, but also enabled the extraotion of the 
underlying physioal strategies. This suooess eneourages the pursuit of other quantum research 
games, as well as dynamic games in other fields. Our future work will focus on both extending 
the gamifioation strategy and also the more general classifioation methodology to new types of 
oontrol problems sueh as those presented in Refs. d^ . d25l and d^ . For all these different 
physioal interaetions we expeot duration-oonstrained problems to exhibit oomplex optimization 
landsoapes. Here players are expeeted to provide a global overview beyond typioal loeal opti¬ 
mization. Another interesting topio of future researoh will be the exploration of how the players 
form strategies. Analyses of the player data will enable us to identify important features of 
quantum optimal eontrol problems and allow an efficient dimensionality reduction. This will 
form the first steps in the efficient training of modern maehine learning algorithms. As a first 
step in this direction, we have collaborated with cognitive researehes to create a new game 
ealled Quantum Minds also available at www.scieneeathome.org. 
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Figures 


(a) 



Figure 1 | The BringHomeWater (BHW) challenge as seen hy the player. The atom is 
represented by the square of its wave funetion, \^jJ{x,t)\‘^, shown as a green liquid. The blue 
eurve represents the potential felt by the atom. Panel h The eontrollable tweezer is initially 
on the left and the atom is trapped in the right statie potential. The player eontrols the optieal 
tweezer by moving a eomputer eursor, pieks up the atom (h), and drags it baek to the target area 
(c), marked by a eyan reetangle, to eolleet points. 
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(a) 




Duration - T 


Figure 2 | Fidelities of the transport problem for different solution durations. Panel a 
shows a subset of the solutions found by the players as turquoise dots. The curves show the 
best solutions found by computer optimizations for each duration. The optimizations shown 
are KASS (red), shovelling (blue), tunneling (yellow), and HILO (purple). The diamonds mark 
the shortest duration with optimal fidelity {F > 0.999) for each optimization method. Panel 
b shows the sweeps from seeds that are generated by players (yellow and blue) and computers 
(red). Player solutions divide into shovelling (blue) and tunneling (yellow) clans. 
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Figure 3 | Shovelling and tunneling clans. The shovelling (tunneling) elan is marked blue 
(yellow) throughout. Panel (a) is a distanee map, showing the distanees between CHOP solu¬ 
tions as defined by Eq. O' Boundaries of the elans are marked with eoloured squares. Panel 
b shows the average trajeetories followed by the elans as thick lines, and a single standard de¬ 
viations thereof as translucent areas. Trajectories are divided into the tweezer amplitude (top) 
and tweezer position (bottom). Panel b shows the wave function fraction (population) in the 
different instantaneous energy eigenstates. 
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(a) T > Tqsl (b) T < Tqsl 



Figure 4 | Optimization landscapes. Panels a and b show the 2D rendering of the high¬ 
dimensional optimization landseape for proeess durations T = 0.40 and T = 0.17, respeetively. 
Green areas mark the spaee probed by CHOP solutions. Panel b is the low-dimensional HILO 
landseape. 
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Methods 


Quantum Moves game interface 

In the online computer game Quantum Moves, the potential created by the optical tweezer 
is represented along the horizontal (position) and vertical (energy) axes (see Fig. ^). The 
potential is updated as the player uses the computer cursor to control the position anointensity 
of the optical tweezer. The sloshy liquid seen by the player on top of the tweezer potential 
is the probability density of the atom position. The probability density (|'0(x, f)p) is defined 
by the quantum mechanical wave function The wave function is updated using the 

time-dependent Schrodinger equation, according to the player-controlled potential. Thus, the 
time evolution is both observed and affected in real time by the player. The time is scaled 
approximately by a factor of 3 ■ 10"^, so that atomic evolution at microsecond time scale is 
experienced over a typical 10 second execution time of the individual game. 

The aim of the games is to find the fastest possible path that transfers the atom into the (sta¬ 
tionary) ground state with unit fidelity, F = 1, in the specified target region. When the tweezer 
is moved quickly, the atomic probability distribution begins to visibly slosh in a fluid-like fash¬ 
ion and the players may readily adopt strategies to prevent or control excitation (sloshing) of the 
wave. To encourage players to search for fast solutions necessary for realistic quantum com¬ 
putations, and to probe the region around the QSL, we introduce a time penalty in the game 
structure. For a thorough discussion on the Quantum Moves player demographics, including 
their scientific backgrounds, we refer the reader to [ISTl and for a similar analysis of Zooni- 
verse.org players see [l32ll . 


KASS 


We tried different seeding and optimization strategies, and the most successful was a Krotov 
algorithm using sinusoidal seed functions and sweeps over the total duration (KASS). KASS 
applies a sweep to an initial random seed Ui{t), where Ui{t) = A{t) is the tweezer amplitude 
and U 2 {t) = x{t) is the tweezer position. The initial random seed is 



( 2 ) 


n=0 


where N is the number of time slices in the path and 6t = 0.002 is the time discretization. Wi is 
the motion at constant speed to the initial position of the atom and back again for A{t) = —100. 
The amplitudes X[n] were selected with a random sign and a norm decreasing as a function of 
the frequency. We generated roughly 2400 seeds using the random summation series and apply 
the Krotov algorithm to these. The next step in KASS is to apply a sweep to the random seeds 

Ui{t)\ 

1. The optimized solution for the tweezer depth and position is linearly contracted in time 

by 5t, A{t),x{t) —)■ A{a ■ t),x{a ■ f) for 0 < a < 1. 

2. The Krotov algorithm is applied using the contracted solution as seed. 

3. Repeat 1-2 until a minimal duration of T = 0.07 is reached. 

Each of the 2400 optimizations took 6 hours of calculation individually. 
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Schrodinger equation in dimensionless units 


The time-dependent Sehrodinger equation solved in Quantum Moves has the dimensionless 
form, 


2 


+ V{x,t)'ip 



where the position x is measured in units of = A/2 = 532nm. This is a typieal period 
of optieal lattice potentials used in experiments with atoms Il20ll . The potential energy 
V{x, t) is given in units of = 1.70peV, about 1/5 of the so-called recoil energy in the lattice 
potential. Time is given in units of E = h/Eu = 0.39ms. In these dimensionless units, our 
optical tweezer is parametrized as 




tweezer 


= A exp I 


2.0(x — Xq)" 

Wo 


where Xq is the position and wl is the waist. In the BHW challenge we have chosen Al = 130 for 
the depth of the static tweezer potential and wo = 0.25 for the width of both tweezer potentials. 
The time-dependent Schrodinger equation is solved using the split-step method. 


Identification of solution clans 

The distance between solutions was calculated using integral of Eq. Q for T ~ 0.17, the 
duration for which the two best CHOP families in Fig. intersect. The solutions were then 
sorted by picking one state at random, and choosing the next state on the list to be the closest 
one in the sense of Eq. O- This process is iterated with the n -f 1-st solution chosen as 
the one closest to the n-th solution. The reachability plot of Extended Data Fig. shows the 
distances between the n -f 1-st and nth solutions. Valleys in the reachability plot clearly identify 
blocks of closely spaced solutions, constituting our solution clans. Clans (valleys) were selected 
by setting a minimum clan size to 200 and an upper threshold 0.05 for the distance between 
consecutive solutions in a clan. These clans are the clans marked in yellow and blue in Fig. 
1^ and Extended Data Fig. As the initial solution is chosen at random, this gives N distinct 
ways of sorting the solutions. We found that the large clans were essentially independent of the 
choice of the initial solution used for the sorting. 


Physical interpretation of solution strategies 

In the tunneling clan, marked yellow in Figs. and and Extended Data FigJ^ the atom is 
collected by tunneling it into the tweezer potential, wich is placed on the lennand side of 
the static potential. As illustrated in Fig. [^, around the time t = 0.06, all of the ~ 500 
yellow solutions move to a very particular location in space, which we interpret as the position 
maximizing the tunneling rate. Instead, the depth of the tweezer potential at this position is only 
weakly correlated with the final fidelity of the solution (Fig[^), permitting large variations in its 
values. This strategy fails for short durations (T < 0.22) as mere is no longer time for the atom 
to tunnel completely between the potential wells. Although seemingly simple, this strategy is at 
variance with the intuition obtained from similar tunneling-based problems [I33ll34l[^ . which 
require careful resonant matching of initial and final state energy levels. In contrast, the CHOP 
solution utilizes a deep potential for the transport tweezer. Combined with precise motional 
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control of the tweezer, this loads the atom direetly into a state with high energy spread (see Fig. 

In the shovelling elan, marked blue in Figs. |^and ^ the tweezer is moved past the position 
of the atom. Therefore the overlap of the tweezer and trie statie potential forms a strong potential 
gradient that aeeelerates the atom towards the target (see Fig. [^). This strategy transfers the 
atom into a superposition of many different instantaneous energy levels with a large energy 
spread AE (see Fig. |^) allowing fast motion in Hilbert spaee towards the target state. As 
illustrated in the inset of Extended Data Fig. |^, this shovelling strategy stays optimal for 
shorter durations than the tunneling strategy. Bemw T = 0.20, however, it eannot be exeeuted 
due to a physieal bound imposed on the speed of the tweezer, and the yellow tunneling elan 
beeomes superior. 


Mapping high-dimensional landscape to 2D 

The visualization of the eontrol landseape maps all solutions to a 2D surfaee. The mapping 
aims to represent distanees Dj^i between any two solutions. This distanee is a Manhattan type 
distanee, 

. JO 

where are resealed tweezer position and amplitude to intervals of unit length. The solutions 
are mapped to the 2D-landseape by the eonstruetion rules: 

1. A randomly ehosen solution is plaeed at the origin {xi, yi) = (0, 0) 

2. The two solutions elosest to the initial solution are given Euelidean eoordinates {xi,yi) 
sueh that Euelidean distanees between them mateh the Manhattan type distanees Di^ 2 , 
Di 3 , and T) 2 , 3 - The points define a triangle with an arbitrary orientation. 

3. The solution with the smallest Manhattan distanee to the previous ones is given eoor¬ 
dinates {xk,yk) sueh that the 2D distanees are as elose as possible to the values Di^k 
i = 1,2,k — 1. This is done with the Nelder-Mead algorithm using the eost funetion. 


Sk= 


(3) 


where denotes the Euelidean distanee between the eoordinates of two solutions in 
2D landseape. 

A 3D figure is obtained by plotting {xi,yi,Fi) where F* is the fidelity of the i-th solution. 
This proeedure gives a qualitative impression of the optimization landseape. Eandseapes are 
eonstrueted from player seeds, random eoiMuter seeds and optimized solutions at the durations 
T ~ 0.17 and T = 0.40 are shown in Eig. Q The multitude of spikes appears beeause gradient 
based optimization leads to nearly vertieaHines in the landseape. This highlights the failure 
of loeal optimization algorithms to explore extended parts of the global landseape. A similar 
landseape for HIEO is also shown in Eig. 


14 




HILO 


HILO, or Heuristically Initialized Local Optimizer, applies the Krotov algorithm to a seed found 
by using the same heuristic as the best CHOP solutions. The motion of the tweezer in a tun¬ 
neling (yellow) and shoveling (blue) CHOP solution can be reasonably approximated by three 
movements at constant horizontal speed and rate of change of the potential depth (see Extended 
Data Fig. J^), which connects three points Pi = {xi, Ai,ti), i = 1,2,3. These three lines 
correspondto the first quick movement to a point Pi, from Pi a slow backwards motion to 
P 2 and a final move to P 3 inside the target area. This defines an optimization problem with 
dimension D = 9. P 3 was set equal to the last point in te,. We noticed that the duration of 
the quick movement (fi) should be as small as possible, reducing the seed space dimension to 
D = 5. This dimension could be further reduced to P = 3 by only changing the position of the 
tweezer using straight lines and letting the amplitude decay exponentially to ^ = —150 after 
ti- Optimizations showed that A = —150 was the optimal value. To reach the desired state 
with high fidelity snaking the tweezer is crucial. The KA is excellent at introducing the snaking 
to stabilize the wave function so this parametrization deliberately does not include these. We 
generated seeds for KA from this space using a direct search. 

The Direct Hilbert Velocity 

In numerous studies the fidelity below the QSL has been found to follow a universal sin^(|T/Tqt 
behavior [[ 8 l. This section briefly summarizes how an expression for the fidelity can be obtained 
using the concept of direct Hilbert velocity Q 0. State |'0(O)) is evolved in time under a con¬ 
trolled Hamiltonian 'H(f) towards a target state \x{T)) at time T. The Krotov algorithm exploits 
the fact that the target state can be propagated backwards in time using the (adjoint) time evo¬ 
lution operator. Using |x(i ^))7 the process fidelity can be evaluated at any instant of time. 


F = \{x{TmT))\^ = \{xmitW- 


It is useful to define the component |0 of the backward evolved target state |x(f))> which is 
orthogonal to the instantaneous state |' 0 (t)). 



The norm of this component reduces at a rate 


Q{t) = ReU{t)\'ip{t) 


Im(^(f) P(f) 


which is denoted the direct Hilbert velocity Q. As described in dH, the trade-off between 
process duration and achievable fidelity obeys the relation 



(4) 


where {Q)t is the time average of Q over the process duration. Equation Q is always true for 
a uniform extension of time and it is valid for any extension of time for an optimal process [|9l- 
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Integrating equation Q for an optimal elan with vanishing initial fidelity leads to 

F = sm‘^(^£{Q)T'dry (5) 

Implying that whenever {Q)t is independent of time, Tqsl = 'k/2{Q)t and for shorter dura¬ 
tions, F{T) = sin^(|T/TQSL)- However, a priori there is no reason that the average direet 
Hilbert speed should be eonstant over an extended range of proeess durations. {Q)t has been 
ealeulated for BHW for the optimal CHOP paths using the method deseribed above, shown 
in Extended Data Fig. [^. Note that |.^) and henee {Q)t is only defined for durations where 
F < 1. Extended DataTigure]^ shows that {Q)t is not eonstant in time, whieh explains the 
deviations from sin^(|T/TQSL) behaviour in Extended Data Fig. |^. How the fidelity drops 
when approaehing the QSE for the different strategies ean be seen in the insert in Extended 
Data Fig. [^. 
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Extended Data Figures 



Duration - T 

Extended Data figure 1 | Deviations from sin^ - behaviour. Colouring as in Fig. 1. Panel 
a shows the fidelity as a funetion of duration for the optimal families resealed by the appar¬ 
ent QSL. The Tqsl is found for eaeh solution by fitting sin^ (aT + b). The blaek line shows 
sin^(7r/2T) for referenee. Panel b shows the direet Hilbert speed {Q)t for the best solutions 
found by HILO, tunneling and shoveling elans respeetively. Note that {Q)t is only defined for 
durations with F < 1 so the curves end at different durations. The varying direct Hilbert speed 
explains the deviation from sin^(|T/T qsl) on the a. 
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Extended Data figure 2 | The reachability plot for Fig. |^. The distance between subsequent 
solutions as calculated by Eq. Valleys identify closely spaced solutions, denoted clans. 
The valleys corresponding to tunneling and shoveling clans are marked with yellow and blue 
respectively. The black line marks the threshold for the distance between consecutive solutions 
in a clan at 0.05. 



18 

























































Extended Data figure 3 | Solutions from CHOP and HILO. The amplitude and the position 
of the tweezer as a function of time for the best player solutions in the tunneling (yellow) and 
shoveling (blue) clans and the best HILO (purple). The dashed purple line shows the initial 
seed used by the best HILO solution (note that the dashed and solid purple lines in the top panel 
overlap). The total duration is T = 0.15. 
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